Coherent light scattering from a two-dimensional Mott insulator 
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We experimentally demonstrate coherent light scattering from an atomic Mott insulator in a 
two-dimensional lattice. The far- field diffraction pattern of small clouds of a few hundred atoms 
was imaged while simultaneously laser cooling the atoms with the probe beams. We describe the 
position of the diffraction peaks and the scaling of the peak parameters by a simple analytic model. 
In contrast to Bragg scattering, scattering from a single plane yields diffraction peaks for any 
incidence angle. We demonstrate the feasibility of detecting spin correlations via light scattering by 
artificially creating a one-dimensional antiferromagnetic order as a density wave and observing the 
appearance of additional diffraction peaks. 

PACS numbers: 42.50.Ar, 67.85.-d, 42. 25. Fx, 37.10. Jk 



Ultracold atoms in optical lattices have become a use- 
ful tool to simulate static phases and the dynamical 
responses of quantum many-body systems [1 j. Recent 
interest has focused on reaching sufficiently low tem- 
peratures and entropies to observe magnetically ordered 
quantum phases [2]. In this context, light scattering has 
been proposed as a new tool to detect these quantum cor- 
relations. Spin correlations could be mapped onto cor- 
relations of scattered light [3 or be detected via diffrac- 
tion peaks from the additional scattering planes for spin- 
dependent probe light [4 . Light scattering would allow 
to measure the temperature of fermions in an optical lat- 
tice [5] or the density fluctuations across the superfluid- 
to-Mott-insulator transition [6j [7] . Since the amount of 
scattered light is usually very small, several proposals 
involve a cavity for the detection [8 j. Without cavities, 
elastic Bragg scattering has been used to demonstrate the 
long range periodic order of thermal atoms in an optical 
lattice despite very low filling factors j9j [10] . It allowed 
the measurement of a change of the lattice constant from 
the backaction of the atoms [9j [10] , their localization dy- 
namics [11] and temperature [10]. Bragg scattering was 
also studied in a far-detuned one-dimensional lattice [12] . 

Here we show coherent light scattering from an atomic 
Mott insulator (MI) in a two-dimensional square lattice 
structure. Scattering from a 2D geometry differs signifi- 
cantly from usual Bragg scattering, because the momen- 
tum transfer needs to be an integer multiple of a recipro- 
cal lattice vector only within the plane. Therefore diffrac- 
tion peaks appear for any angle of the incoming beam, 
which is experimentally more convenient. In our setup, 
we use five probe beams in a molasses configuration that 
simultaneously laser cool the atoms. Each of these mo- 
lasses beams yields distinct diffraction peaks in the far- 
field images. We quantitatively compared the diffraction 
patterns with model calculations and confirmed the co- 
herent nature of the scattering process. We artificially 
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FIG. 1. Schematics of light scattering, (a) Light diffraction 
from atoms in a ID lattice. An incoming wave with wave vec- 
tor ki is diffracted under an angle 0. (b) Resulting differential 
scattering cross section (#), as given by Eq. (pi for our ex- 
perimental parameters. Shown are the cases for unity filling 
of the lattice and atom number N = 16 (blue curve) and for 
a z-Neel antiferromagnet along one direction (N — 28), from 
which only one spin component is detected (green curve) . The 
grey shaded area indicates the opening angle of our imaging 
system, (c) Experimental setup. Atoms in a 2D optical lat- 
tice are illuminated with four in-plane molasses beams. In 
situ and far-field diffraction images are recorded with a high- 
resolution optical imaging system. 



prepared ID antiferromagnetic order as a density wave 
and observed additional diffraction peaks, thus demon- 
strating the usability of light scattering for the detection 
of global spin correlations. 

We begin by introducing a simple analytic ID model 
to illustrate the underlying physics. Atoms on a ID 
lattice (lattice spacing a\ at = Ai at /2, Ai at is the lattice 
wavelength) are driven by an incoming light field (wave- 
length A) from the x direction [Fig.[TJa)] with wave vec- 
tor lq = ke x , where k = 2tt/X and is the unit vec- 
tor along x. The scattered light is detected at a point 
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r, defined by the angle 6 with the z-axis, such that 
r = r e r = r sin e x + r cos e z . The position of the l-th 
atom is = I a\ at e x and its distance r\ to the detection 
point is in far-field approximation T\ = |r— xj| w r— x/-e r . 
In our model, each atom emits a spherical wave, which 
at the detection point can be written as 
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Here, / denotes the coherently scattered field amplitude, 
Si =k.i • xj is the phase imprinted by the incoming light 
field, and K = k — with the wave vector k = k e r 
in the observed direction. The differential cross section 
^f(K) oc | e" lK ' Xi | 2 is obtained by summing over the 
field amplitudes from all N x atoms. As a result, we ob- 
tain the angular dependence of the scattered light field, 
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with distinct maxima when the field amplitudes of neigh- 
boring atoms interfere constructively, i.e. K-(x/— xj+i) = 
2tt • n, where n is an integer that denotes the diffraction 
order. The height of the diffraction peak is proportional 
to N% whereas the peak width scales as 1/N X . The angles 
n , under which the diffraction maxima can be observed, 
are given by 



sin n = 1 
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The trivial case n = gives the forward scattered light 
(#o = 90°) > independent of a\ at and A. For our ex- 
perimental parameters (A = 780 nm, ai at = 532 nm), 
Eq. @ can be additionally fulfilled only for n = — 1, 
yielding the corresponding minus first diffraction order 
at 0_i = -27.8° and 207.8°. These two out of plane 
scattered waves ensure the momentum conservation in 
the z direction. Fig.[TJb) (blue curve) shows a polar plot 
of j@(0), displaying the forward scattered light and the 
two peaks, one of which is captured by our imaging sys- 
tem (gray shaded region). If only every second lattice 
site is occupied (e.g. after removing one spin component 
in an antiferromagnetically ordered sample), the period- 
icity of the system is doubled. In this case, there are two 
possible diffraction orders (in the upper half plane) at 
0AFM = 15 5 o and afm = 6> _ i = _27.8° [green curve in 

In our experiment, we prepared 2D Mis of 87 Rb atoms 
in an optical lattice. The atoms were confined in a single 
antinode of a vertical standing wave. Two pairs of hori- 
zontal laser beams provided the square lattice structure 
(see Ref. [13]). In order to detect the atoms in the opti- 
cal lattice, the lattice depth was increased to ~ 300 /iK, 
thereby freezing the atom distribution. The atoms where 
then illuminated with an optical molasses, which con- 
sisted of two pairs of retroreflected laser beams oriented 
along the horizontal lattice axes [see Fig.Qc)]. A fifth 



molasses beam entered from the reverse direction of the 
imaging system (not shown in Fig.[TJc)). The molasses 
was red detuned with respect to the free space resonance 
by 45 MHz and the total scattering rate was 60 kHz. We 
detected the fluorescence photons with a high resolution 
microscope objective with half opening angle a = 43°. 
The objective can be moved in the z direction by a dis- 
tance Az < 100 /im between the atom position and the 
focal plane within 50 ms using a piezo scanning device. 

Our experimentally obtained diffraction images are 
shown in Fig.[2|a) for four different distances Az and an 
illumination time of 200 ms for each image. For Az = 
we observed the in situ atom number distribution, con- 
sisting in this case of 147 atoms in a MI shell with unit 
occupancy and diameter of 6 /im. For larger Az, we ob- 
serve the build up of the far-field distribution with dis- 
tinct diffraction peaks. We compared the experimental 
data with a numerical calculation of ^| (K) using the ac- 
tual atom distribution of the image at Az = [Fig.|2jb)]. 
For this purpose, we coherently summed over all spher- 
ical waves Fi(ri) emitted by the atoms with phases Si 
given by the incident driving fields. Our model assumes 
that all four horizontal molasses beams are diffracted in- 
dependently. A spherical wave for the emission pattern 
is used, because the different local polarizations in the 
molasses result in all possible orientations of the atomic 
dipole. The calculated far-field distribution is in good 
qualitative agreement with the experimental images. Our 
simulation only includes the coherently scattered light, 
whereas the experimental data also shows a significant 
incoherent background. 

For a more quantitative analysis, we recorded diffrac- 
tion patterns of Mis for different atom numbers [Fig.[3ja)- 
(c)]. We evaluated cuts (angular sectors of width 4°) 
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FIG. 2. Light scattering from a 2D ML (a) Experimental 
images from the same atomic sample for four different dis- 
tances Az of the focal plane with respect to the atom position. 
The lowest image shows the in situ atom number distribution, 
whereas the upper image shows the far-field diffraction pat- 
tern, (b) Simulated diffraction patterns obtained from a 2D 
numerical model (see text for details) , using the reconstructed 
atom distribution (lower image) from (a). Red arrows indi- 
cate the directions of the optical molasses beams. 
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FIG. 3. Analysis of the diffraction patterns, (a)-(c) In situ 
images (N =126, 199 and 279 atoms) and the correspond- 
ing far-field diffraction pattern (Az = 40 /im) from the same 
experimental run, with illumination time 500 ms each, (d) 
Angular distribution of the differential scattering cross sec- 
tion obtained from cuts in (c) through the diffraction peaks 
(dark blue) and the background signal (light blue). The peak 
is fitted with a Gaussian (red line), (e) Peak height A versus 
N together with a quadratic fit. (f) Peak width w versus N, 
with a fit to w oc 1/ ^/N . (g) Total background signal <t& (light 
blue) as obtained from the constant background in (d) and 
total signal in the peaks a p (dark blue, scaled by a factor of 
10) obtained from the fits to the diffraction peaks (see text). 

through the diffraction peaks and the background sig- 
nal [see Fig.|3jc)] and applied an appropriate coordinate 
transformation in order to obtain the angular distribu- 
tion || (0), as shown in Fig.knd). We fitted the resulting 
peaks with a Gaussian (height A, 1/y/e width w, center 
position and offset fixed at the background value). 
The peak position, averaged over all experimental runs 
is = 27.4(6)°, in excellent agreement with the ex- 

pected value of \0-i\ = 27.8°. The error is dominated 
by the systematic uncertainty of ±1 /im in the determi- 
nation of Az. The peak height scales quadratically with 
the atom number [Fig.[3je)], illustrating the coherent na- 
ture of the diffraction peak. The peak width scales as 



w oc 1/y/N [Fig.|3[f)l, in agreement with the result from 
the ID model [Eq. ij2|], assuming N x = y/~N atoms in 
one dimension. 

There are several mechanisms that lead to a back- 
ground outside the diffraction peaks. The first mecha- 
nism is the deviation from integer occupation of the lat- 
tice sites caused by density fluctuations in the system 
[6]. The second mechanism is the spread of the atoms 
in their potential wells, which reduces the peak height 
via the Debye- Waller factor. In our case, the dominant 
mechanism, however, is inelastic light scattering in the 
molasses configuration. We used the far-field images to 
extract the power scattered into the detected diffraction 
peaks and the power scattered into the background. In 
the cuts outside the diffraction peaks, we found a con- 
stant background (^f )b and calculated the total incoher- 
ent scattering cross section — 4^(^)5, assuming an 
isotropic intensity distribution. The total scattering cross 
section from the five molasses beams can be estimated as 
a p « 10 sin(6Li)27rAiu 2 cos(0_i). The factor cos ac- 
counts for the ellipticity of the diffraction peaks due to 
the effective ellipticity of the atomic cloud, when viewed 
under the angle 6-\. Fig.[3|g) shows and <r p , which 
both scale linearly with the number of atoms. From the 
slopes, we find a fraction f p = a p / (<r& + a p ) « 3% of the 
power scattered into the detected peaks. Only about 20% 
of the coherently scattered light is diffracted into these 
peaks, while the forward scattered part is not detected 
here. 

In addition to the four diffraction peaks from the hori- 
zontal molasses beams, a fifth weaker peak is clearly visi- 
ble in the center left part of the far-field images [see white 
circle in Fig.[3|c)]. This peak results from the diffraction 
of the molasses beam which is shone in from the direction 
of the imaging system. It shows that our single plane of 
a few hundred atoms in the optical lattice acts as a "mir- 
ror" for the incoming laser beam. 

Finally, we demonstrated that light scattering can be 
used for the detection of spin correlations. As an exam- 
ple, we created a ID z-Neel antiferromagnetic order along 
the x direction of the lattice using our recently demon- 
strated single-site addressing technique [14] . We sequen- 
tially flipped all atomic spins in every second row of the 
lattice from F = 1 to F = 2 before we applied a resonant 
laser that removed all atoms in F = 2. Figs. [4^ a) and (b) 
show the resulting fluorescence image in the focal plane 
together with the reconstructed atom number distribu- 
tion. The corresponding experimental and theoretical 
diffraction patterns are displayed in Figs.pFc) and (d). 
The two predicted diffraction peaks of -l s ^and -2 nd or- 
der along x are clearly visible in the experimental picture, 
although our atomic sample consisted of only 57 atoms. 
We obtain the usual peak position of = 27.7(6)° 

along the y direction, whereas the two peaks along x 
are found at \6 X _ X \ = 14.5(6)° and \6 X _ 2 \ = 27.4(6)° [see 
Figs.|2fe),(f)], in good agreement with the expected val- 
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FIG. 4. Light scattering for a ID antiferromagnetic order 
in the density, (a) The atoms in every second row were re- 
moved, (b) Reconstructed atom number distribution from 
(a), (c) Resulting far- field image (Az = 25 /xm) with two 
diffraction orders in the x direction, (d) Simulated diffraction 
pattern using the atom number distribution of (b). (e),(f) 
Angular distribution of the differential scattering cross sec- 
tion obtained from cuts along the x and y directions together 
with Gaussian fits. 

lies. We prepared a ID antiferromagnetic order, because 
the additional diffraction peaks that would arise for a 2D 
antiferromagnetic order lie outside the opening angle of 
our imaging system. However, the position of the diffrac- 
tion peaks could be varied in a 2D geometry by changing 
the angle of the incident beams. An alternative is to use 
shorter wavelength probe light, e.g. near-resonant with 
the 5S - 6P transition at 420 nm for 87 Rb. 

Our results could be extended to the study of vari- 
ous density [6[ [7] or spin [3j |4] correlations in optical 
lattices. Most proposals suggest weak non-destructive 
probing, which restricts the signal to only a few pho- 
tons per atom. In our alternative approach, we projected 
the correlations onto the density before the detection by 
means of an optical molasses, which yields a signal of 
thousands of photons per atom. This is possible for den- 
sity correlations, as e.g. the number squeezing in a MI 
can be mapped on the mean density by parity projection 
P~3j[15]. For spin correlations, we have demonstrated the 
feasibility of removing one spin state and observing the 
additional diffraction peaks from the density structure. 
This avoids spin selective coupling of the probe light to 
the atoms [4], which is incompatible with simultaneous 
laser cooling. In the 3D case, antiferromagnetic order al- 



lows scattering from an additional plane, but it requires 
careful alignment of the probe beam angle to match the 
Bragg condition. In contrast, light scattering from a 2D 
system yields additional peaks from the same incident 
beam, which is also convenient for the extraction of the 
spin correlation length from the relative height or width 
of the diffraction peaks. Finally, we note that the detec- 
tion of spin correlations via light scattering does not rely 
on the high aperture of our imaging system. The diffrac- 
tion peaks could also be detected in a restricted angular 
range. 
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